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Semiannual  Technical  Report  on 
Atmospheric  Pressure  Gas  Lasers 

covering  the  period 
June  1,  1973-November  30,  1974 


SUMMARY 


I.  Amplification  of  Two  High- Intensity  Nanosecond  TEA.  Laser  Pulses 


Yongyut  Manichaikul  studied  the  effect  of  vibrational  relaxation 
upon  the  replenishment  of  the  upper  laser  level  when  two  high  inten¬ 
sity  pulses  are  passed  through  a  CC^  TEA  amplifying  medium.  He 
showed  that  with  increasing  pressure  in  the  amplifying  medium,  the 
large  signal  gain  of  the  second  pulse  increases  until  it  reaches  a 
plateau  corresponding  to  partial  recovery  of  the  laser  inversion  due 
to  redistribution  of  the  vibrational  population  associated  with  the 
asymmetric  stretching  mode.  The  experiment  indicated  that  the  001 
level  of  CO2  was  repopulated  by  the  higher  Onm  levels  at  the  rate 
that  was  directly  proportional  to  the  CO2  partial  pressure. 

Details  of  the  experiment  are  described  in  Appendix  I,  which  is 


a  reprint  of  the  RLE  Quarterly  Progress  Report  No.  Ill  of  October 
15,  1973. 


II.  Mode  Locking 

C.  P.  Ausschnitt  has  begun  experimental  work  aimed  at  a  better 
understanding  of  active  and  passive  mode  locking  of  CO2  lasers. 
Concurrently,  work  is  in  progress  to  obtain  closed- form  expressions 
for  mode-locked  pulses  produced  both  actively  and  passively.  The 
literature  on  the  subject  is  very  large;  yet  in  spite  of  the  vast 
amount  of  work  on  the  subject,  there  still  is  no  compact  simple 
theory  that  would  both  predict  the  steady- state  mode-locking  solu¬ 
tions,  and  also  be  amenable  to  the  study  of  the  stability  of  such 
steady- state  solutions,  their  transient  behavior,  their  susceptibility 
to  noise,  and  their  stability. 

Considerable  progress  has  been  made  in  an  effort  to  accomplish 
this  program.  A  recent  paper  by  Kurokawa  on  injection  locking  pro¬ 
vided  the  basis  for  an  analysis  which  is  both  simpler  and  more 
flexible  than  analyses  presented  heretofore.  A  paper  on  the  theory 
of  forced  mode  locking  forms  Appendix  II.  While  the  steady- state 
solutions  worked  out  in  this  paper  have  been  obtained  previously  in 
the  literature,  we  present  here,  for  the  first  time,  the  solution 
of  the  forced- mode- locking  problem  of  an  inhomogeneous ly  broadened 
medium  and  study  the  transient  behavior  of  the  buildup  of  mode- locked 
pulses,  as  well  as  the  stability  of  the  various  mode-locked  "super 


nu.  3es" ,  both  for  the  homogeneously  and  inhomogeneously  broadened 
laser  medium.  Work  is  currently  in  progress  that  utilizes  the  same 
formalism  to  analyze  saturable  absorber  mode  locking. 


III.  Closed-form  Analysis  of  Electronic 
Distribution  and  Pumping  Rates 


The  work  carried  out  on  this  topic  has  appeared  in  the  February  1974 
issue  of  the  Journal  of  Applied  Physics;  pp.  781-791. 
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Reprinted  from 

Quarterly  Progress  Report  No.  Ill,  October  15,  1973 


(VI.  APPLIED  PL  ASMA  RESEARCH) 


3.  P.  K.  Cheo  and  R.  L.  Abrams,  "Rotational  Relaxation  Rate  of  C02  Laser  Levels," 
Appl.  Phys.  Letters  1_4,  47-49  (1969). 

4.  K.  F.  Herzfeld  and  T.  A.  Litovitz,  Absorption  and  Dispersion  of  UP ra sonic  Waves 
(Academic  Press,  Inc.,  New  York,  1959). 

2.  AMPLIFICATION  OF  TWO  HIGH-INTENSITY  NANOSECOND 
TEA  COz  LASER  PULSES  (AHN) 

National  Science  Foundation  (Grant  GK-37979X) 

U.  S.  Army  -  Research  Office  -  Durham  (Contract  DAHC04-72-C-0044) 

Y.  Manichaikul 
Experiment 

We  have  previously  reported  on  the  generation  and  amplification  of  high -intensity 
nanosecond  pulses.1  Two  or  three  of  these  pulses  were  produced.  They  wore  from 
the  P(  1 6)  transition,  2  ns  wide  (FWHM),  separated  by  12  ns.  When  these  pulses  were 
focused  into  a  three-electrode  laser  amplifier  as  shown  in  Fig.  VI-i2,  a  peak  intensity 
of  2-3  MW/cm2  was  obtained.  A  beam  splitter  was  used  so  that  the  intensity  of  the 
pulses  could  be  monitored.  The  input  and  output  detectors  were  as  shown  in  Fig.  Vl-12 
In  this  experiment  the  detected  input  signals  were  delayed  100  ns  by  using  60  ft  of 


R  100*  FLAT  MIRROR 


FROM  OSCILLATOR  OF  HIGH  -  INTENSITY  m  PULSE 
GENERATING  SYSTEM  DESCRIBED  PREY,  DUSLY.1 


Fig.  VI-12.  Experimental  arrangement  for  amplification  of 
high-intensity  ns  pulses.  (See  Y.  Manichaikul.1) 

RG-8  cable.  The  add  mode  of  a  Tektronix  oscilloscope  was  used  to  display  the  signals 
for  both  input  and  output  pulses  on  the  same  screen.  The  two  detectors  were  calibrated 
against  each  other  by  comparing  the  oscilloscope  picture  of  the  input  and  output  pulses 
without  discharge  exciting  the  three -electrode  laser  amplifier.  Figure  \I-1  3a  shows 
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(b) 

Fig.  VI-13.  Output  and  input  of  AHN  expe  iment.  Total  pressure: 

200  Torr.  Gas  mixture:  CO.:N_:He  =  X:4:100. 

^  ^  2 

Intensity  of  input  pulse:  0.75  (MW/cm  )/div.  Inten- 

2 

sity  of  output  pulse:  1.12  (MW/cm  )/div.  Time: 

5t  ns/div.  (a)  Amplifier  off.  (b)  Amplifier  on. 

the  oscilloscope  display  for  this  case. 

In  order  to  probe  the  temporal  evolution  ot  the  three -electrode  laser  amplifier,  we 
first  fired  the  amplifier  and,  after  a  chosen  delay  time,  the  oscillator.  In  general,  the 
oscillator  was  fired  ~30  ps  after  the  onset  of  the  discharge  in  the  amplifier  for  two  rea¬ 
sons:  first,  we  wished  to  avoid  the  effects  on  our  measurements  of  the  shock  waves 
generated  by  the  discharge.  Second,  we  wished  to  be  certain  that  the  symmetric 

stretching  (SS)  and  bending  (B)  modes  of  CO-  had  equilibrated  with  each  other  at 

^  2 

slightly  above  the  kinetic  temperature  of  the  gas. 

Measurements  on  the  amplification  of  high -intensity  ns  pulses  were  made  at  200  Torr 
of  COz:N2:He  mixtures.  The  ratio  of  these  mixtures  was  C02:N2:He  =  X:4:100,  where 
CO,  partial  pressure  was  varied  from  3.  5  to  35  Torr  partial  pressure.  Small-signal 
gain  of  this  three -electrode  laser  amplifier  in  each  case  was  measured  by  a  cw  CO, 
laser. 

Results 

Figure  VI-1 3b  illustrates  the  input  and  output  pulses  vhen  the  amplifier  is  turned 
on.  Four  such  measurements  were  made  and  their  average  was  taken  at  each  CO, 
partial  pressure  studied.  We  have  found  that  the  RG-8  cable  used  for  the  time  delay 
introduces  some  distortion  in  the  input  signals.  This  distortion  can  be  accounted  for 
if  the  first  (second)  pulses  of  the  input  and  output  pulses  from  the  amplification  mea¬ 
surements  are  compared  with  the  first  (second)  pulses  of  the  input  and  output  pulses 
when  the  amplifier  was  evacuated. 


INPUT 

OUTPUT  / 


(a) 


INPUT 

OUTPUT 


(LSG,/I 


Fig.  VI-14.  LSG2/LSGj  vs  partial  CO. 
pressure. 
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Table  VI-1.  Experimental  results. 


No. 

— 

pc°2 

(Torr) 

SSG 

±10% 

LSGj 

±10% 

Ta 

vv 

(°  K) 

AN 

u 

N  (0) 
u 

1 

3.  5 

0.  50 

0.  30 

1180 

0.  12 

2 

6.  5 

0.  70 

0.  36 

1150 

0.  11 

3 

11.0 

1.  50 

0.  76 

1155 

0.  14 

4 

20.  0 

1.00* 

0.  52 

845 

0.  11 

5 

23.  0 

1.  50* 

0.75 

860 

0.  13 

6 

24.0 

2.  00* 

0.96 

890 

0.  14 

7 

27.  0 

2.  50* 

1  1 5 

935 

0.  15 

8 

35.  0 

2.  00 

0.  81 

830 

0.  10 

Notes:  Not  measured  directly;  calculated  from  LSGj  and  the 

peak  intensity  of  the  pulse. 


bcG  =  — —  is  the  small-signal  gain  across  the  tube, 

in  2 

Here  the  intensity  is  less  than  1  W/cm  . 

LSGj ,  large -signal  gain  of  the  first  pulse. 

Ta  ,  temperature  of  the  asymmetric  stretching  mode  cal- 
vv  culated  from  SSG. 

AN 

- — ,  fractional  depletion  of  the  00°  1  population  by  an 

Nu(0)  ns  pulse. 

AN  ,  calculated  from  the  large -signal  gain  and  the  inten- 
U  sity  of  the  pulse. 

N  (0),  obtained  from  SSG. 
u 
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Figure  VI-14  shows  LSG^LSGj  s  the  partial  pressure  of  C02  studied.  We  have 

,  A  *outi  *ini 
-  - r - , 

Ain  i 

where  i  =  1,2,  with  1  and  2  representing  first  and  second  pulses.  LSGj  (LSG2)  is 
the  large-signal  gain  of  the  first  (second)  pulses.  The  following  observations  can  be 
made  from  these  measurements,  (i)  LSG2/LSGj  is  less  than  unity.  This  is  to  be 
expected,  since  the  first  pulse  had  depleted  a  fraction  of  the  population  from  the  00°  1 
level  of  C02.  (ii)  The  ratio  LSG^LSGj  is  approximately  0.  8  at  C02  partial  pressure 
of  3.  5  Torr  and  the  ratio  increases  slowly  to  0. 9  as  CC>2  partial  pressure  increases 

to  20  Torr  or  higher,  which  is  as  expected,  since  the  00°  1  level  of  CCL  was  being 

jp  3  ^ 

repopulated  by  the  higher  O  n  m  levels  at  a  rate  that  is  directly  proportional  to  the 

C02  partial  pressure. 

Table  VI-1  gives  other  experimental  results  of  interest.  We  found  that  the  large- 

signal  gain  of  our  pulses  is  approximately  one -half  the  small-signal  gain,  and  the  frac- 

AN 

tional  depletion  of  the  00°  1  level,  - —  ,is  between  0.  10  and  0.  15. 

N  (0) 

u 

A  theoretical  model  for  the  amplification  of  high-intensity  nanosecond  pulses  is  being 
developed.  We  shall  present  the  theory,  and  make  a  comparison  of  theory  and  experiment 
in  a  future  report. 
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APPENDIX  II 


A  THEORY  OF  FORCED  MODE  LOCKING* 

** 

by  H.  A.  Haus 


Abstract: 

A  new  theory  of  forced  mode  locking  is  deve2.oped  in  the 
frequency  domain.  Replacement  of  the  discrete  cavity  mode 
spectrum  by  a  continuum  reduces  the  problem  to  the  solution 
of  an  eigenvalue  problem  of  a  total  differential  equation.  We 
treat  the  mode  locking  of  both  a  homogeneous  and  an  inhomo¬ 
geneous  laser  medium.  The  equations  of  mode  locking  by  modu¬ 
lation  of  resistance  and  reactance  in  a  cavity  containing  a 
homogeneously  broadened  laser  medium  are  those  of  the  harmonic 
oscillator  of  quantum  mechanics.  The  "supermodes"  of  the  mode 
locked  system  as  functions  of  the  continuum  frequency  are  the 
eigenfunctions  of  the  harmonic  oscillator.  A  stability  analysis 
of  this  system  shows  that  only  the  lowest  order  supermode  is 
stable  with  respect  to  small  perturbations.  The  properties  of 
this  supermode  are  the  same  as  those  previously  derived  by 
Kuizenga  and  Siegman.  The  same  analysis  lends  itself  to  the 
treatment  of  mode  locking  of  a  cavity  containing  a.i  inhomoge- 

*  Work  supported  by  Joint  Services  Electronics  Program  (Contract 
DAAB07-71-C-0300) .  US  Army  Research  Office  -  Durham  (Contract 
DAHC04-72-C-0044) . 

**  Electrical  Engineering  Department  and  Research  Laboratory  of 
Electronics,  Massachusetts  Institute  of  Technology,  Cambridge, 
Mass.  02139. 


neously  broadened  laser  medium.  Again  one  finds  that  the 
lowest  order  supermode  with  a  single  peak  (as  a  function  of 
frequency)  is  stable,  all  higher  order  modes  are  unstable. 
The  transient  build-up  of  a  supermode  is  analyzed.  Effects 
of  modulator  detuning  are  treated. 


V/l/l 


Introduction 


The  literature  on  mode  locking  is  large,  comprising  both 
theoretical  and  experimental  investigations! 1] .  The  theoretical 
investigations  may  be  roughly  divided  into  two  groups,  (a)  those 
which  treat  the  mode  locking  process  in  the  time  domain,  e.g. 
the  excellent  paper  by  Kuizenga  and  Siegman[2] ,  and  (b)  those 
which  treat  the  phenomenon  in  the  frequency  domain  (notably  the 
paper  by  McDuff  and  Harris).  [3]  In  this  paper  we  present 
a  useful  extension  of  the  mode  locking  analysis  in  the  frequency 
domain.  When  setting  up  the  analysis  of  saturable  absorber  mode 
locking  by  a  somewhat  novel  approach  we  soon  discovered  that  this 
approach  lent  itself  very  well  to  a  better  understanding  of  the 
forced  mode  locking  phenomenon  and  permitted  a  broader  investi¬ 
gation  of  the  influence  on  mode  locking  of  changes  in  various 
parameters. 

This  paper  is  devoted  to  forced  mode  locking.  In  partic¬ 
ular  the  following  issues  are  addressed:  1.  Mode  locking  of  q 
homogeneously  broadened  laser  via  single  side-band-pair  genera¬ 
tion,  using  both  amplitude  and  phase  modulation.  2.  Tuning  of 
the  mode  locking  element  drive  away  from  the  frequency  corre¬ 
sponding  to  the  cavity  mode  spacing.  3.  Changes  in  position  of 
the  mode  locking  element  along  the  cavity.  4.  The  difference 
of  mode  locking  between  lasers  possessing  inhomo -reneous  versus 
homogeneous  broadening.  5.  Transient  bu^ld-up  of  a  mode  locked 


(2) 


pulse.  6.  Stability  of  the  mode  locked  pulses. 

We  consider  the  problem  in  the  frequency  domain  utilizing 
the  theory  of  injection  locking  as  recently  described  by  Kurokawa. [4] 
Mode  locking  differs  from  injection  locking  only  in  that  the  in¬ 
jected  signal  is  generated  internally  to  the  oscillator  cavity  by 
the  signal  present  within  the  cavity.  A  complication  of  the  laser 
mode  locking  is  that,  if  it  is  to  be  of  any  use,  many  cavity  modes 
must  be  locked  and  hence  the  analysis  calls  for  the  solution  of 
many  simultaneous  equations.  This  fact,  however,  can  be  bent  to 
one's  advantage.  As  long  as  the  modes  under  consideration  are 
numerous,  one  may  treat  them  as  a  continuum  and  replace  the  dif¬ 
ference  equation  by  a  differential  equation.  This  approximation 
has  the  advantage  that  the  eigenvalue  techniques  of  differential 
equations  can  be  brought  to  bear  on  the  mode  locking  problem. 

The  approximation  of  replacing  a  discrete  set  of  modes  by  a 
continuum  implies  that  the  pulses  predicted  by  this  theory  must 
be  separated  (in  time)  sufficiently  so  as  not  to  overlap.  As 
soon  as  overlap  of  the  pulses  becomes  important,  the  replacement 
of  the  discrete  spectrum  by  a  continuum  is  not  legitimate. 
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II.  Mode  Locking  as  Injection  Locking 

! 

4 

Injection  locking  of  a  cavity  with  a  single  mode  containing 
a  negative  resistance  may  be  analyzed  with  the  aid  of  the  equiv¬ 
alent  circuit  of  Fig.  1.  Here  the  impedance  “Z  (A)  represents 

d 

the  active  impedance  which  is  in  general  complex  and  a  function 

of  the  amplitude  of  the  complex  circulating  current  I.  If  a 

i 

voltage  source  E  is  inserted  between  the  two  impedances  as  in- 

dicated,  then  the  equation  for  this  injection  locked  oscillator 

| 

is  given  by 

I 

F  =  tZc(w)  -  Zg(A)]  I  (2.1) 

In  the  sequel  we  shall  follow  Kurokawa's  notation  who  uses  I 
for  the  complex  current  amplitude  and  A  for  its  (real)  magni¬ 
tude. 

I 

Eauation  (2.1)  can  be  adapted  to  the  case  cf  a  mode  locked  cav 
itjy  involving  the  interaction  among  many  axial  modes.  We  comprise 
the  current  amplitudes  of  all  the  modes  of  interest  in  a  current 
amplitude  column  vector  I.  If  the  modes  in  the  cavity  are  un¬ 
coupled,  in  the  absence  of  a  mode  locking  crystal,  then  the  im¬ 
pedance  matrix  of  both  the  passive  and  the  active  part  of  the 
circuit  is  diagonal.  We  shall  simplify  the  analysis  by  assuming 
that  the  cavity  impedance  may  be  described  by  the  diagonal  matrix 

Rc  f  1  +  j  xc(u,i)]  (2.2) 

where  Rc  is  the  equivalent  cavity  series  resistance  assumed 

identical  for  all  modes;  x  contains  the  normalized  reactive 

•**  c 

components  of  the  impedances  of  each  mode,  a  diagonal  matrix; 
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to  the  extent  that  the  reactances  may  differ  from  mode  to  mode, 
x  is  not  proportional  to  the  identity  matrix.  We  indicate  by 
cj .  the  frequency  of  the  injected  signal  (as  well  as  oscillation) 

i 

for  a  general  cavity  mode.  The  frequencies  of  the  driven  oscil¬ 
lations  of  all  other  modes  are  "spaced"  from  the  center  frequency 
by  a  multiple  of  wm,  where  is  the  frequency  of  modulation. 

(See  Fig.  2) 

Mode  locking  produces  sidebands  for  each  of  the  oscillating 
amplitudes.  The  generalization  of  (2.1)  to  this  case  is 

M  I  =  (  1  -  Z  (A)  +  j  x  (w.)  )  I  (2.3) 

z  (A)  is  the  normalized  impedance  matrix  contributed  by  the 

A,U 

active  medium,  a  diagonal  matrix,  but  not  proportional  to  the 
identity  matrix.  On  the  left  hand  side  appears  the  matrix  M 

++'•4 

which  is  a  nondiagonal  matrix  and  accounts  for  the  generation  of 
sidebands  by  mode  locking.  The  detailed  form  of  the  matrix  M 
depends  upon  the  nature  of  the  mode  locking  scheme. 

It  should  be  pointed  out  that,  in  accordance  with  Kurokawa, 
we  use  series  equivalent  circuits  to  represent  the  cavity  modes 
interacting  with  the  laser  medium.  If  E-field  amplitudes  of 
the  laser  modes  were  identified  with  voltages  and  H-field  am¬ 
plitudes  with  currents,  a  parallel  equivalent  circuit  would  result. 
The  reader  has  then  the  choice  of  either  (a)  adhering  to  this 
identification  and  reading  "admittance"  instead  of  "impedance" 
everywhere  in  the  sequel,  or  to  identify  currents  with  E-field 
amplitudes  and  voltages  with  H-field  amplitudes. 


(5) 
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III.  Modulation  of  Resistance. 

Suppose  the  cavity  contains  a  resistive  medium  which  is 
modulated  in  time,  so  that  its  resistivity  has  the  time  depen¬ 
dence: 


r* 


h 


constant  +  (  1  -  cos  u>  t) 

m 

If  the  modulated  medium  is  positioned  inside  the  cavity  so  that 
all  modes  have  the  same  field  distribution  through  the  medium 
(i.e.  the  mode  locking  element  is  relatively  thin  and  is  near 
one  of  the  mirrors)  then  the  medium  contributes  to  the  resis¬ 
tance  R  of  the  cavity  a  time  dependent  component 
c 


R  (  1  -  cos  tii  t) 
m  m 


v'here  the  remaining  time  independent  contributions  may  be  lumped 
wi th  R  . 

°  jwkt 

A  current  |  IjJ  e  pertaining  to  the  k-th  mode,  flowing 

through  P.  ,  produces  a  voltage 
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The  (equivalent  injection  locking)  voltage  at  the  freauency 
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of  the  k-th  mode  is  therefore: 
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(3,2) 


The  modulation  matrix  of  (2.3)  becomes,  in  the  present  case: 
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Introducing  (3.3)  into  (2.3)  one  obtains  a  second  order  dif¬ 
ference  eouation. 

We  may  now  make  a  aeries  of  approximations  so  as  to  enable 
us  to  find  simple  solutions  to  this  equation.  First  of  all,  if 

t 

many  modes  are  to  be  locked,  and  if  the  variation  of  the  gain 
over  a  frequency  range  corresponding  to  the  mode  spacing  is 
small,  then  the  difference  equation  can  be  replaced  by  a  second 
order  differential  equation.  One  has 
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where  I,  the  column  vector  containing  the  (complex)  current 
amplitudes  of  the  modes,  is  now  replaced  by  a  continuous  function 


of  freouency  deviation  from  material  line-center, 

is  the  frequency  separation  of  the  spectral  lines. 

now  identified  by  the  continous  variable  to  *=  lim  / 

k  •*■<» 

k  is  the  mode  number  and  assumes  both 
positive  and  negative  values.  Defining  the  modulat 
R 

M  «  -52—  one  obtains  for  (2.3)  [5]: 
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The  modes  are 

>  to  ,  where 
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:ion  coefficient 


(3.5) 


In  the  sequel  we  specialize  to  specific  cases. 


IV.  Homogeneously  Broadened  Laser  Linewidth. 


When  the  laser  is  homogeneously  broadened,  then  the  impe¬ 
dance  contributed  by  the  active  medium  z  ,  is  proportional 
to  1/(1  +  P/P  )  where  P  is  the  total  power  in  the  laser 

S 

cavity  in  all  modes  and  Pg  is  the  saturation  power.  This  is 
true  over  the  entire  linewidth,  i.e.  for  all  values  "  of 
the  function  z  (A,  w) .  We  shall  make  \-he  further  approxima- 
tion  that  the  gain  profile  (linesnape)  of  the  laser  medium,  and 
hence  the  frequency  dependence  of  negative  resistance  r,  can 

"  a 

be  replaced  by  a  parabola.  This  assumption  holds  provided  the 
loss  line  is  reasonably  close  to  the  peak  of  the  gain  curve  so 
that  no  modes  oscillate  within  the  frequency  range  in  which  the 
line  shape  of  the  medium  differs  appreciably  from  the  parabolic 
approximation.  With  these  assumptions  we  have 


z  (A,  w)  =  r 
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(4.1) 


where  wM2  is  related  to  the  curvature  at  the  peak  of  the  gain 
line  and  xs  a  measure  of  the  laser  medium  bandwidth.  The  de¬ 
pendence  on  power  P  of  r  is 
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(4,2) 


where  rQ  is  the  small  signal  value  of  r. 


We  have  added  a  reactive  contribation  to  the  later  medium  which , 
by  the  Kramers-Kronig  relations,  must  be  associated  with  a  fre¬ 
quency  dependent  real  part. 

Introducing  (s.l)  into  (3.5)  we  obtain 
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r  /  1  -  —  \  +  rj  —  +  j  x  (u) 
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(4.3) 


In  (4.2)  the  assumption  was  made  that  the  gain  of  the  medium 
does  not  get  modulated  by  the  passing  laser  pulses;  the  medium 
gain  is  pulled  down  on  a  time  average  basis  only.  This  assump¬ 
tion  amounts  to  the  statement  that  the  population  recovery  time 
(T^)  of  the  medium  is  long  so  that  the  gain  change  is  caused  by 
the  depleting  action  of  many  successive  pulses,  each  pulse  changing 
the  population  slightly.  Gain  modulation  does,  in  fact,  tend  to 
counteract  mode  locking  and  will  be  treated  ir  another  publica¬ 
tion.  Hence  the  theory  developed  in  chis  paper,  like  the  theory 
of  Siegman  and  Kuizenga,  is  good  only  for  relatively  long  T^. 

Next  consider  the  meaning  of  x  (w) .  This  is  the  cavity 

C 

mode  reactance  as  a  function  of  mode  number,  here  indicated  by 
the  continous  variable  w.  x_  is  zero  if  the  mode  is  excited 

C 

at  resonance  and  is  ncnzerc  for  a  frequency  deviation  6w  from 
mode  resonance  and  given  by 

X  .  2  —  Q 
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4 


(4.4) 


(10) 


where  C  is  the  quality  factor  of  the  cavity  mode  and  wQ  is 
the  resonance  frequency  (which  may  be  considered  to  be  a  constant 
*or  all  modes  within  ':he  narrow  laser  line).  Eq.  (4.3)  is  par¬ 
ticularly  simple,  if  we  require  that 


xc  (w)  +  r  —  =  0 

This  implies  that 

6w  „  rw  1 
“o  “M  2(3 


(4.5) 


namely  that  the  frequency  deviation  6w  from  mode  resonance  is 

proportional  to  the  frequency  separation  from  the  material  line 

center.  Such  a  condition  can  be  achieved  by  proper  adjustment 

of  the  frequency  of  the  modulator.  Indeed,  we  may  write 

w  *  k  w_  where  n  is  the  order  of  the  mode  as  measured  from  ma¬ 
in 

teriai  line  center.  Therefore,  if  the  modulator  frequency  wra 
is  set  so  that 


6w  *  k(w  -  Aw) 
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we  find  for  w  : 
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(4.6) 


In  this  way,  we  have  achieved  the  desired  frequency  deviation 
as  a  function  of  mode  number  k;  (4.6)  implies  that  the 
modulation  frequency  is  set  equal  to  the  mode  resonance  frequency 
separation  as  modified  by  the  laser  medium.  With  (4.5)  satis¬ 
fied,  Eq.  (4.3)  has  a  pure  real  coefficient  on  the  right  hand 
side.  The  resulting  equation  is  identical  with  the  harmonic 
oscillator  equation  of  quantum  mechanics,  where  I  plays  the 
role  of  the  wave  function  and  w  the  role  of  the  spacial  vari¬ 
able. 

We  can  now  bring  to  bear  the  entire  formalism  of  the 
quantum  mechanical  harmonic  oscillator  on  the  problem  of  loss 
mode  locking  in  a  homoaoneously  broadened  medium.  The  solutions 
of  (4.3)  are  the  well  known  Hermite  polynomials  [6] 

(O  =  Hn(a  exp  -I  c2 

where  the  variable  F,  is  a  spacial  variable  normalized  in  terms 
of  the  parameters  of  the  harmonic  oscillator. 


The  role  of  the  harmonic  oscillator  quantity  -fS2/2m  is  played 
by  the  parameter  Mw^2  and  the  sprinq  constant  has  to  be  iden¬ 
tified  with  2r/wM2.  With  this  identification  we  have 
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where 


n  = 


"[ M  ~  /m  oT“ 

\  /  -  v  m  M 


(4.8) 


The  eigenvalues  of  energy  of  the  harmonic  oscillator  follow 
the  well  known  law  hv  (n  +  1/2)  where  n  is  an  integer.  The 
role  of  energy  is  played  by  the  parameter  r  -  1  in  equation 
4.3.  We  find  for  the  quantization  of  this  parameter 


r  .  i  B  2  /Mr  %  (n  +  1/2) 


(4.9) 


Equation  (4.7)  defines  " supermode s"  or  "mode  locking  modes" 
of  the  mode  locked  oscillator.  The  higher  the  order  n  of  the 
"supermode \  the  more  structure  in  the  frequency  spectrum.  Since 
the  Fourier  transform  of  the  spectrum  (4.7)  leads  to  the  same 
time  dependence,  we  find  that  mode  locked  solutions  may  exist  of 

higher  and  higher  temporal  structure  exhibiting  more  and  more 
pulses.  (See  Fig.  3) 

For  a  given  excess  gain,  rQ  >  1,  and  given  modulation  M, 
the  bandwidth  of  the  lowest  order  supermode  (n  «=  0)  is  pro¬ 
portional  to  ,  i.e.  the  square  root  cf  the  pressure  in  a 

homogeneously  broadened  gas  laser.  Equation  (4.9)  determines 
the  power  in  the  mode  locked  cavity  for  the  various  "supermodes". 
Indeed  by  solving  for  r  from  (4.9)  and  then  using  (4.2)  we 
find  the  following  relationship  for  the  power 


Pg  1  +  V  +  /2y  +  v 


(4.10) 
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where  y 


-  (2n  +  l)2 
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In  the  limit  of  a  small  modulation  coefficient  M,  (4.10) 
becomes 

P  w 

r  =  ro  -  1  -  2  ro^  ~  <"  +  1/2)  (4.10a) 

Ps  “M 

We  see  that,  for  a  given  modulation  coefficient,  the  power  de¬ 
creases  with  "supermode"  order.  Mother  interesting  fact  may 
be  noted  from  (4.9).  The  highest  order  supermode  supportable 
in  a  given  cavity  of  specified  small  signal  gain,  is  obtained 
by  setting  r  =  rQ  in  (4.9).  We  note  that  the  smaller  M,  the 
modulation  coefficient,  the  higher  the  allowed  order.  This 
appears  somewhat  counter-intuitive  until  one  realizes  that  we 
have  not  discussed  which  supermodes  are  stable  and  thus  WDuld, 
in  fact,  exist  in  the  steady  state.  We  have  merely  found  that 
weak  modulation  leads  to  weak  selectivity  among  the  supermodes 
and  hence  to  presumably  a  noisy  output.  The  larger  the  modu¬ 
lation  coefficient  M,  the  fewer  supermodes  are  allowed.  The 

question  of  stability  of  the  suDermodes  will  be  taken  up  in 
Section  XI. 

The  bandwidth  is  given  roughly  by  the  turning  points  of 
the  Hermite  solutions,  namely  by  2u>e  where 
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or,  using  (4.8)  and  (4.10) 
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(4.11) 


The  bandwidtn  is  proportional  to  the  4-th  root  of  the  modulation 
coefficient  and  the  square  root  of  n  +  1/2.  Our  analysis  is  ap¬ 
plicable,  as  long  as  (2tog/tow  )  >>  1,  or 


2  \l  2  /M7F  ^  (n  +  1/2)  »  1 
which  is  not  too  severe  a  restriction  as  long  as  wM/w  »  1. 


(4.12) 
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V.  Phase  Modulation 


Next,  we  investigate  mode  locking  via  modulation  of  a 
reactance.  Suppose  the  reactance  X^  is  modulated  so  that  it 
has  a  time  dependence 


X  [  1  -  cos  w  t] 
m  m 


(5.1) 


If  one  defines  the  modulation  coefficient 


M  = 


(5.2) 


one  obtains  an  equation  identical  to  (4.2)  except  that  now  M 
on  the  left  hand  side  of  (4.2)  has  to  be  multiplied  by-j. 
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d2 1 
dw2 
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f1  “  -4V rj  — +  -xc(w) 


(5.3) 


This  equation  can  be  reduced  to  the  form  of  (4.2)  through  a 
change  of  variables.  One  defines  the  new  independent  variable 


w"2  =  u2//~j 


(5.4) 


Kith  this  change  of  variable,  one  finds 


I 


(5.5) 


Clearly  the  aolutiona  of  this  equation  are  the  same  as  before. 
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j)fi  exp - -  (1  -  j)ft2 
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(5.6) 


where  the  normalized  frequency  ft  has  been  defined  in  (4.8). 
The  phase  of  the  mode  locked  complex  current  amplitudes  varies 
across  the  frequency  profile.  Using  the  fact  that  the  eigen¬ 
values  of  (5.5)  are  the  same  as  the  ones  previously  obtained  in 


(4. 9)  one  has 


(r  -  1)  (1  -  j)  -  (1  *  j)  L  +  r  fi-j  .  2  ^  % (n  +  V2)  ? 
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Equating  the  imaginary  parts  of  this  equation,  one  finds  the 
amount  of  the  detuning  of  the  cavity  modes  from  resonance. 


(r-  1)  »  -  |xc+r~ 


We  find  that  the  modulator  frequency  wm  obeys  a  slightly 
different  condition  than  (4.6),  if  solution  (5.6)  is  to  apply. 
Indeed,  using  (4.4)  we  find 


x  =2  —  Q  =  2 - 6w  +  k  (cu  -  Aw)  =  -r  — -  -  (r  -  1) 


(5.8) 


where  6cj q  is  the  frequency  deviation  of  the  oscillation  from  the 
resonance  frequency  of  the  mode  at  laser  medium  line  center  (k  =  0) 
To  the  extent  that  Aw  -  ^  /  0,  this  frequency  deviation  changes 
as  one  goes  away  from  line  center.  Matching  the  k-dependent  parts 
of  (5.8) ,  one  has: 
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which  is  identical  to  (4.6).  Hence,  solution  (5.6)  requires  that 
the  modulation  frequency  be  tuned  to  the  mode  separation  frequency, 
as  modified  by  the  reactance  of  the  laser  medium.  Matching  the* 
k- independent  parts  of  (5.8)  one  has: 


(r  -  1) 


The  oscillations  occur  off  mode  resonance  (as  modified  by  the 
laser  medium).  Equating  the  real  parts  of  (5,7),  one  finds 
a  relation  for  the  power  of  oscillation. 


n/2/mT  %(n  +  1/2) 

S  -  -  —  -  -  -  - 


We  have  found  that  phase  modulation  works  similarly  to 
amplitude  modulation,  but  that  the  cavity  modes  do  not  oscil¬ 
late  at  cavity  line-center. 
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A  comparison  of  the  pulses  produced  by  amplitude  modula¬ 
tion,  or  phase  modulation*  can  be  made  if  one  studies  the  in¬ 
verse  Fourier  transform  of  the  spectra.  We  concentrate  solely 
on  the  lowest  order  supermode.  Taking  into  account  that 
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we  may  make  the  following  identifications: 
a)  AM  modulation 


M 

F“m  “m 


b  =  0 


so  that  the  exponential  multiplied  by  the  carrier  becomes 


exp  -  f!  M 7r  wMt2)  exp  j«0t 


The  pulsewidth,  defined  as  the  width  between  the  half  power 
points  is 


x  =  2/a  ln5  =  2 


\l  ln2  - i 
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h)  FM  modulation 


The  exponential  is 

exp  -  -  /M/2r  t2  exp  j(wQt  -  -  /v/2r  u>m  t2] 

2  2 

The  width  of  the  phase  mode  locked  pulse  is  by  -v  2  greater 
than  that  of  the  amplitude  mode  locked  pulse.  This  is  in  agree¬ 
ment  with  a  result  due  to  Kuizenga  and  Siegman  obtained  by  a 
different  method,  if  their  notation  is  brought  into  correspondence 

with  our  notation.  .  , 

In  addition,  the  frequency  of  the  phase-modulation  pulse 

is  chirped,  having  the  "instantaneous  frequency 
-  /m72t  urn  wM  t 

The  frequency  decreases  with  time  tor  the  assumed  modulation 


of  reactance  (5.1). 
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VI.  Modulator  Detuning,  Resistance  Modulation 


We  shall  now  suppose  that  the  modulator  frequency  u)m  is 
not  set  equal  to  the  frequency  spacing  of  the  cavity  modes  (as 
determined  in  the  presence  of  the  laser  medium).  Therefore, 
the  modes  at  greater  "distance"  from  line  center,  become  pro¬ 
gressively  more  detuned.  The  reactance  xc  of  a  cavity  mode 
at  a  deviation  6<u  from  a  center  frequency  wQ  is  given  by 
(4.4).  If  w  is  not  set  equal  to  the  separation 
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(6.1) 


of  the  cavity  resonance  frequencies  in  the  presence  of  the  laser 
medium,  the  function  in  the  brackets  on  the  right  hand  side  of 
(4.3),  i.e.  the  impedance  of  modes  and  laser  medium  as  a  function 
of  mode  number  k,  becomes  complex .  One  has 


x  +  r  —  «*  —  k(u)  -  Aw)  +  r  — - 

c  m 

“M  Uo  wM 
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If  one  defines  the  detuning  parameter  6  by 


where  Awc  is  the  cavity  mode  bandwidth  c  .o  obtains: 
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Introducing  this  dependence  into  (4.3),  we  find  the  differential 
equation 
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Again  it  is  not  difficult  to  find  a  solution  for  this  differential 
equation  using  a  shift  in  origin  of  the  independent  variable. 
Introducing  the  shifted  frequency  to'  by  the  relation 
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we  may  obtain  the  new  differential  equation 


(6.6) 


Again,  this  is  the  harmonic  oscillator  equation.  Note,  however, 
that  the  shift  of  origin  of  the  independent  variable  is  imaginary. 
The  lowest  order  supermode  can  be  shown  to  be  time  delayed  (with 
respect  to  the  modulator  phase)  by  the  amount: 
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Writing  down  the  eigenvalue  relation  for  this  equation  one  obtains 


62wm2  2 /Mr  to  (n  +  1/2) 

— -  r=  _ X. _ _ 
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One  finds  that  mode  locking  may  be  prevented  by  excessive  de¬ 
tuning  of  the  modulation  freouency.  The  detuning  introduces 
a  term  analogous  to  a  loss  increase  of  the  cavity. 

The  effect  of  detuninq  may  be  interpreted  as  follows: 
Minimum  loss  of  pulse  would  be  achieved  if  the  pulse  passed  the 
mode  locking  crystal  at  the  time  of  its  peak  transmission.  The 
repetition  rate,  however,  is  inconsistent  with  the  cavity  reso¬ 
nance  freouency  spacinq .  Hence,  the  cavity  reactances  are  non¬ 
zero.  This  calls  for  a  reactive  contribution  of  the  mode  lockinq 
crystal  which  is  achieved  when  the  peak  of  the  pulse  passes  at 
a  time  different  from  the  peak  transmission  time.  Additional 
loss  due  to  the  modelocking  crystal  is  then  unavoidable. 
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VII.  Changes  in  Position  of  Mode  Locking  Crystal 


So  far  we  have  assumed  that  the  crystal  was  near  one  of 
the  mirrors  so  that  all  cavity  modes  had  roughly  the  same  field 
through  toe  crystal.  The  equivalent  circuit  for  a  cavity  con¬ 
taining  a  resistivity  p  and  having  normalized  magnetic  field 
patterns  of  modes  k  and  l,  Kk  and  respectively,  pos¬ 

sesses  a  coupling  resistance  Rk£  between  the  modes  k  and 
£  of  the  form 


/  P  h. 


ht  dv 
hk dv 


(7.1) 


where  the  integrals  are  taken  over  the  cavity  volume. 


Now  consider  an  equivalent  current  Ik  exp(jwkt)  of  the 

k-th  mode  and  ask  for  the  voltage  at  the  freouency  k  +  1.  This 

is  given  by  the  couplina  resistance  R  which  has  the 

K  t  if  K 

followina  dependence  upon  the  field  patterns 


jwkt  /  p  (t)  h  •  h,_  dv  jw,t 
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Suppose  that  we  have  standing  waves  inside  the  cavity  and  all 
modes  of  interest  have  roughly  the  same  cross-sectional  depen- 


a 


dence.  In  this  case,  one  finds  for  the  integral  in  the  limit 
of  a  crystal  of  short  length  L,  positioned  at  a  distance  X( 
from  one  of  the  mirrors. 


x  +  L 
o 


lim 

Aw  L/c  •+•  0  x 


j  dx  cos  8k  x  cos  8k  +  ^  x 


Rk+1,  k  x  +  L 


(7.3) 


J  dx  cos2  +  ^  x 


A  to 

COS  -  X 


where  Bk  =  —  =  («  +  kAw)/c 

Xk 


We  thus  find  that  the  strongest  coupling  is  achieved  when  the 
position  of  the  crystal  is  near  one  of  the  mirrors  and  no  coupling 
is  achieved  if  the  crystal  is  positioned  at  the  center  of  the 
cavity.  This  finding  depends  in  part  on  the  assumption  that 
the  mode  locking  crystal  produces  only  a  single  pair  of  sidebands. 
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VIII.  Higher  Order  Sideband  Generation 

The  generalization  to  include  the  case  in  which  more  than 
a  s  ngle  pair  of  sidebands  is  generated  is  not  difficult;  what 
is  difficult  is  the  solution  of  the  resulting  equation.  Indeed 
if,  say,  two  pairs  of  sidebands  are  generated  by  the  modulation, 
located  symmetrically  around  the  modulated  carrier,  then  the 
second  derivative  of  the  current  amplitude  has  to  be  replaced 
by  a  superposition  of  a  fourth  derivative  and  a  second  derivative. 
The  resulting  equation  is  much  more  difficult  than  the  one  obtained 
thus  far.  Whereas  solutions  for  equations  of  this  form  have  been 
obtained  in  the  literature [ 7] ,  they  are  complicated  and  have  not 
yet  been  explored  with  the  present  purpose  in  mind. 

Another  way  of  generalizing  the  equations  is  to  replace  M  — - 

dw2 

by  an  operator.  An  integral  operator  is  found  to  be  appro¬ 
priate  for  the  analysis  of  a  saturable  absorber  and  will  be  taken 
up  elsewhere. 
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IX.  Inhomogeneous ly  Broadened  Medium 


If  the  medium  is  inhomogeneous ly  broadened ,  then  the  power 
in  any  particular  cavity  mode  determines  the  pulldown  of  the 
gain  line  at  that  particular  cavity  frequency  (provided,  of  course, 
that  the  separation  of  the  cavity  modes  is  large  compared  with 
the  homogeneous  linewidth  of  the  medium) .  The  adaptation  of 
(4.3)  to  the  case  of  an  inhomogeneously  broadened  medium  is  ac¬ 
complished  by  changing  the  chrract  r  of  the  dependence  of  the 
negative  resistance  upon  intensity .  For  simplicity  we  omit  the 
reactive  contribution  of  the  laser  medium,  an  assumption  justi¬ 
fied  by  the  fact  that  the  overall  gain  profile  tends  to  get 
"flattened"  by  the  simultaneous  oscillations  of  many  cavity  modes 
and  the  individual  modes  oscillate  at  the  center  of  the  local 
hole.  The  eouation  is: 


I! 


d2  I 
dw2 


■•r*  „ 

1  -  -  +  3  x  (w) 


1  + 


V 


(9.1) 


where  A  in  the  denominator  is  the  square  of  the  amplitude  at 
the  particular  frequency  of  interest  and  Ag2  is  proportional 
to  the  saturation  intensity.  (9.1)  is  a  nonlinear  differential 
equation  for  the  current  I  (remember  that  A  is  the  magnitude 
of  the  current).  It  is  a  difficult  equation  to  solve;  yet  it  is 
possible  to  gain  considerable  insight  into  the  nature  of  its  so¬ 
lutions  by  giving  it  a  physical  interpretation.  Let  us  separate 
the  complex  current  amplitude  into  an  amplitude  and  a  phase  factor 


I 


(9.2) 


=  A  e 


j<f> 


Taking  the  second  derivative  of  (9.2)  one  obtains 

ill  =  A  e^  +  2jJ  A  e^  +  j*  A  -  {2  A  e^ 
dw2 


Introducing  (9  2)  and  (9.3)  into  (9.1)  one  has 


*»  » 

M  w  2  (A  -  <f>  2  A) 
m 


1 


1 


A 


(9.3) 


(9.4) 


M  u  2  (aV  +  2Af)  =  x  A  (9.5) 

in  c 

Here  we  have  written  (*)  for  d/dto(  )  in  order  to  emphasize 

f 

the  analogy  with  the  equations  of  motion  of  a  particle  in  polar 
coordinates;  the  frequency  variable  plays  the  role  of  "time",  the 
amplitude  A  is  analogous  to  "radius".  The  force  field  is 
'radius"  and  "time"  dependent.  If  one  assumes  that  all  cavity 
modes  oscillate  at  line  center  of  the  respective  cavity  resonances, 
then  xc  =  0  everywhere  and  the  force  field  becomes  a  central  one. 
In  this  case,  one  may  define  a  potential  function,  the  derivative 
of  which  give;  the  force.  The  potential  function  plotted  against 
"radius"  A  as  a  function  of  "time"  w  is  shown  in  Fig.  4. 

Consider  first  the  case  in  the  absence  of  modulation.  In 
the  particle  motion  analog,  this  corresponds  to  a  massless  parti¬ 
cle.  The  particle  will  always  look  for  a  position  in  which  it  is 
exposed  to  zero  force  field.  As  we  can  see  from  Fig.  4.  ,  the 


potential  hill  travels  as  a  function  of  time  towards  the  oriqin 
and  then  the  origin  becomes  a  hill  in  its  own  right.  If  one 
looks  for  a  node  locked  "supermode",  one  looks  for  the  motion 
of  a  particle  that  starts  with  zero  velocity  at  a  time  t  =  0 
(zero  slope  at  center  frequency) ,  and  then  approaches  the  origin 
(A  =  0)  where  it  comes  to  a  full  stop.  The  equation  of  motion 
for  the  massless  phrticle  is  then  very  simple,  it  stays  at  the 
top  of  the  hill  starting  infinitesimally  to  the  left  and  then 
searches  for  the  position  in  which  the  force  is  infinitesimally 
small  yet  directed  from  right  to  left.  VThen  the  potential  hill 
has  made  it  to  the  origin,  the  particle  has  followed  it  there 
and  then  stays  there  for  ever.  This  case  can  be  solved  very 
simply  analytically  by  setting  '’.he  righthand  side  of  equation 
(9.4)  =  0. 

Consider  next  the  case  when  there  is  modulation  and  when 
the  particle  has-  to  be  assigned  a  nonzero  mass.  Let  us  consider 
at  first  the  case  of  very  small  mass  (modulation).  At  first, 
when  the  force  field  is  strong,  one  may  still  consider  the  par¬ 
ticle  to  seek  out  positions  in  the  potential  field  at  which  it 
is  exposed  to  practically  zero  force.  This  means  that  the  par¬ 
ticle  follows  roucrhly  the  same  motion  as  we  have  determined 
earlier  by  assuming  zero  mass  for  the  particle,  except  when  it 
approaches  the  origin.  Near  the  origin,  the  force  field  ap¬ 
proaches  zero, but  the  particle  has  acquired  a  finite  velocity 
and  kinetic  energy.  In  order  to  come  to  a  full  stop  at  the  origin, 
it  has  to  expend  the  kinetic  energy  by  climbing  up  a  hill.  This 
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means  that,  before  the  particle  comes  to  full  rest,  the  origin 
must  have  already  risen  above  its  surroundings  so  that  a  par¬ 
ticle,  near  the  end  of  its  motion,  will  have  to  climb  up  a  hill. 
Translated  into  the  language  of  the  mode  locked  "mode"  this 
means  that  the  frequency  bandwidth  of  the  mode  is  wider  than  in 
the  absence  of  modulation  (remember  in  the  language  of  the  equiv¬ 
alent  particle  the  origin  of  the  radial  coordinate  develops  into 
a  hill  after  the  massless  particle  has  reached  the  origin) . 
Further,  it  is  clear  that  a  particle  with  mass  requires  an  ini¬ 
tial  force  in  order  to  be  nudged  towards  the  origin,  which  means 
the  particle  has  to  start  slightly  to  the  left  of  the  potential 
hill  at  t  =  0.  Again,  in  the  language  of  mode  locking,  this 
corresponds  to  a  power  of  the  mode  near  line  center  that  is 
smaller  than  the  power  in  the  absence  of  mode  locking.  The  model 
is  useful  also  for  determining  what  happens  when  the  strength  of 
mode  locking  is  increased  excessively.  This  correspond®-  to  a 
very  massive  particle.  The  potential,  as  defined,  does  not  in¬ 
volve  the  mode  locking  amplitude  M  and  hence,  in  the  equiva¬ 
lent  language  of  the  particle,  the  initial  hill  becomes  smaller, 
the  heavier  the  particle.  A  supermode,  symmetric  at  line 
center,  corresponds  to  an  initial  condition  for  the  analog  parti¬ 
cle  starting  at  rest  to  the  left  of  the  potential  hill.  A  par¬ 
ticle  which  starts  initially  at  rest  slightly  to  the  left  of  the 
potential  hill  may  be  too  massive  to  make  it  near  the  origin 
before  the  potential  distribution  has  changed  beyond  recall,  pre- 


venting  the  particle  from  acquiring  sufficient  kinetic  energy. 
The  particle  never  makes  it  to  the  origin.  No  symmetric  super¬ 
mode  exists  for  excessively  high  values  of  M. 

There  are  also  higher  order  steady  state  solutions,  if  the 
modulation  is  not  too  strong.  One  may  start  the  particle  on 
the  down  slope  of  the  hill.  (Note  the  peak  intensity  at  center 
frequency  for  this  mode  is  less  than  for  the  lowest  order  super¬ 
mode)  .  The  particle  may  go  through  the  origin  and  climb  up  the 
hill  on  the  other  side.  After  oscillating  back  and  forth,  the 
particle  can  be  brought  to  rest  at  the  origin.  Such  solutions 
exist  only  for  the  particles  of  finite  mass  (i.e.  finite  modu¬ 


lation).  Equation  (9.4)  has  two  adjustable  aprameters,  rQ  -  1, 
the  excess  small  signal  gain  at  line  center  (w  =  0)  and  the 


normalized  modulation  parameter 


/M  % 
r0  WM 


For  a  range  of  these  two 


parameters,  solutions  are  obtained.  Fig.  5  shows  plots  of  A/A 

s 

vs  normalized  frequency  w/w.,  ,  for  various  choices  of 

/m'  w 

the  parameters  yj—  ~  and  r0  -  1.  Fig.  6  shows  the  range  in 


the  (rQ  -  1), 


-  (0 

M  m 


r0  WM 


plane  over  which  mode  locked  solutions 


are  obtained.  Fig.  7  shows  the  decrease  of  the  oscillation  am¬ 


plitude  at  line  center  A(0)/A  ,  as  a  function  of  normalized 

S 

modulation,  for  different  excess  small  signal  gains. 


X.  Transient  Buildup  of  Mode  Locked  Laser 


So  far  we  have  concentrated  on  obtaining  steady-state  super¬ 
mode  solutions.  The  case  of  the  homogeneously  broadened  laser 
medium  lends  itself  nicely  for  study  of  the  buildup  of  the  mode 
locked  pulses  (supermodes) .  Since  we  have  treated  the  pulse  spec¬ 
trum  as  a  continuum,  and  a  continuous  spectrum  predicts  one  single 
pulse  only,  we  must  ask  how  the  analysis  can  be  adapted  to  the 
description  of  the  buildup  of  a  sequence  of  discrete  pulses. 

We  must  recall  that  the  continuous  spectrum  was  introduced 

in  the  first  place  in  order  to  obtain  an  approximate  solution  to 

the  difference  equation.  In  fact,  the  supermodes  obtained  thus 

far,  functions  of  frequency  w,  must  be  interpreted  as  defining 

a  closely  spaced  discrete  spectrum,  with  spectral  lines  at  kw^, 

with  k  an  integer.  A  perturbation  of  the  frequency  of  each 

axial  mode,  i.e.  replacement  of  kw  by  kw  +  fiw,  ,  implies  a 

m  m  k 

deviation  of  the  oscillation  frequency,  Re  fiw^  ,  and  growth  of  the 
axial  mode  at  the  rate  -Im  fiwj,.  A  transient  analysis  then  entails 
the  introduction  of  the  frequency  perturbation  fiw^  and  an  expan¬ 
sion  of  all  parameters  of  the  fundamental  difference  equation  in 
terms  of  this  frequency  perturbation.  One  may  go  to  the  continuum 
limit,  treating  fiw  as  a  function  of  lim  kw^  w ,  and  carry  out 
the  expansion  on  the  differential  equation  (4.3). 

We  shall  disregard  all  derivatives  of  the  parameters  of  (4.3) 
with  respect  to  fiw,  except  that  of  x  .  This  is  tantamount  to 
disregarding  all  energy  storages  other  than  the  electromagnetic 
energy  storage  of  the  cavity  modes. 

We  now  turn  to  the  expansion  of  x  ,  a  function  of  mode 

c 


\  / 


number  k.  If  the  oscillation  of  an  axial  mode  occurs  at  a  fre- 

ouency  different  than  the  steady  state  frecrurncy,  then  x  (k.  ku  ) 

cm 

has  to  be  replaced  by  x  =  x (k,  k to  +  610,  )  =  x  (ktu  ) 

c  c  in  Jc  o  in 

3xc(k,  oj) 

+  -  or  in  the  continuum  limit 

3(d 


x„  =  x  (w)  +  6o>  x  ' 
c  c  c 


(10.1) 


where  6a>  is  itself  a  function  of  to  koi  and  we  have  set 

m 


x 


c 


f 


3xc(k,  oj) 
3to 


> 

When  the  steady  state  oscillation  frequencies  are  sufficiently  near 

the  resonance  frequencies  of  the  axial  modes,  x  '  can  be  treated 

c 

as  a  constant. 


d2 1 
du2 


+ 


i  x  '  I 
c 


(10.2) 


Clearly  xc'  >  0  by  Foster's  reactance  Theorem.  Equation 
(10.1)  can  be  solved  simply  if  one  assumes  that  the  width  of  the 
medium  / (ruM  ^  ~ ^rowM2 ^ '  an(^  ^ence  the  laser  pulse  width,  is 
independent  of  time.  We  shall  first  assume  that  this  is  the  case 
and  then  show  the  generalization  to  the  case  when  the  time  varia¬ 
tion  of  r  is  taken  into  account. 

Suppose  the  modulation  is  resonant,  w  =  Ato 

'  m 

xc(“)  +  r  —  =  0  in  (10.1).  If  we  set 


P 

s 


1  +  r  —  “)"1 

“m  20/ 


/ 


and  introduce  this  expression  into  (10.2)  we  obtain 


Mw 
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d2I 

dw2 


r  w' 

_o _  -  i  -  0 


1  + 
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M  J 


I  =  j6w  x  *  I 


Settinn 


I  =  exp  -  i  -4  E  ^  %<“> 

,  2  co 


‘/r 


with 


b)  =  VM/r  /w  w 
n  o  m  M 


(10.4) 


(10.5) 


where  TR  is  the  time  interval  between  successive  pulses,  we 

have  /  vq2  (w)  dm  =  1  and  the  enerqy  in  a  sinqle  pulse  is  qiven 

by  A2  T  ,  the  power  P  is  equal  to  A2.  Now  i6w  A  can  be 
K  > 

reqarded  as  the  time  derivative  of  A,  dA/dt,  thinking  of  the  am¬ 
plitude  of  the  Fourier  spectrum  to  be  a  function  of  time. 

Let  us  make  this  assertion  mathematically  precise.  We  have 
been  treat inq  the  mode  spectrum  as  a  continuum.  In  fact,  the 
mode  spectrum  is  discrete,  and  the  continuum  replacement  only 
served  to  obtain  an  approximate  solution  to  the  difference  eoua- 

tion.  Denote  the  Fourier  transform  of  v  (to)  by  v  (t):  Noting 

o  o 

that  we  are  dealing  with  a  periodic  process  of  long  period  TR,  we 
have  for  a  sequence  v  (t)  ,  of  slowly  varying  amplitude  A(t)  : 

T/2 

/  dt  e~jnAwt  A  (t )  v  (t) 
yTw  _fp  o 


1 

/2v 


A(t) 


Tr/2 

/  at  v  (t) 

-V2 


1 

/2? 


00 


A (t )  /  dt  e^wt  v  (t) 

o 

—  00 


1 

/2? 


A (t)  v0  (w) 


The  first  step  is  legitimate,  because  A(t)  is  assumed  to 
vary  negligibly  within  one  period  TR  and  hence  it  can  be  pulled 
out  from  underneath  the  integration.  The  second  step  is  justified 
by  our  assumption  of  negligible  overlap  of  the  pulses.  The  final 
expression  shows  that  we  may  treat  the  Fourier  spectrum  both  as 
time  and  frequency  dependent. 

Returning  now  to  (10.4)  we  obtain,  using  *10. 5) : 


« 

A  = 


i  -  yis; 


A 


(10.6) 


We  may  make  the  van  der  Pol  assumption  (really  a  necessary  as¬ 
sumption  in  view  of  our  disregard  of  the  variation  of  rwuJ)  i.e. 

M 


With  this  assumption,  (10,6)  written  in  terms  of  P  =  A2 
reduces  to 
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o  p2 
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9  g  d 
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where 


(rQ  -  1)  -/Mr 


w  i 

-  m  1 

0  v  i 
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The  solution  to  this  equation  is 


2a  t 

P  e  g 

p  o 

s  1  +  e  q  pq 


(10.9) 


Fiqure  8  shows  a  plot  of  this  solution  vs  normalized  time  2a  t, 

g 

using  the  normalized  power 


a_ 

-3-  P 


We  now  take  into  account  the  change  in  the  bandwidth  of  the 
laser  medium  by  replacing  Equation  (10.3)  with 


i  -  -si| 
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(10.10) 


It  should  be  noted  that  we  assumed  that  the  modulation  fre¬ 
quency  was  equal  to  the  cavity  mode  spacing,  that  the  mode 

locking  was  on  resonance.  To  the  extent  that  w  is  usually 

m  J 

kept  fixed,  but  the  mode  spacinq  varies  with  amplitude  because 
of  the  changing  dispersion  of  the  laser  medium  this  is  an  approx¬ 
imation.  It  can  be  shown,  however,  that  in  most  cases  the  "de- 
tuninq"  of  the  mode  locking  due  to  a  change  of  the  cavity  mode 
spacinq  caused  by  chanqes  in  r  is  negligible.  We  start  with 
the  assumption  that  the  bandwidth  of  the  mode  locked  pulse  is 
a  function  of  time,  in  (10.5)  we  replace  to  2/r  by  w  2/r_ 
where  Wj  is  a  function  of  time.  This  assumed  expression  for 
the  current  wave  form  obeys  the  eigenvalue  equation 


M  to  2 
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d2 1 
dto2 


I  =  - 


(10.11) 


Now  the  time  derivative  of  the  assumed  waveform  is 


dl 

dt 
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(10.12) 

The  time  derivative  of  the  mode  locked  pulse  contains  a  hermite 
polynomial  of  second  order.  This  is  the  consequence  of  the  fact 
that  a  changing  width  (in  frequency  space)  of  the  pulse  requires 


growth  of  the  wings  at  the  expense  of  the  center.  Replacing  the 
right-hand  side  of  (10.2)  bv  the  time  derivative,  utilizing  (10.12) 
and  balancing  the  w  independent  terms  and  the  terms  proportional 
to  to2  one  obtains 
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with  r 


(10.14) 


These  are  two  coupled  first  order  differential  equations  in  the 
variables  Wj  and  A.  Since  they  are  nonlinear,  they  require  in 
general  a  computer  solution. 

We  shall  derive  here  approximate  solutions  to  these  two 
equations.  Multiplication  of  Eg.  10.14  by  A  so  as  to  obtain 
an  equation  for  the  variable  A2  =  P,  subsequent  introduction 
of  the  van  der  Pol  approximation  leads  to 
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(10.15) 


One  may  expand  the  solutions  of  (10.14)  and  (10.15)  in  terms  of 

a  parameter  rf  smallness  /m  to  /w.. .  The  bandwidth  itself,  as  a 

m  M 

function  of  time,  varies  more  slowly  with  time  than  the  power  and 
hence  to  zero  order  one  may  disregard  the  contribution  of  band¬ 
width  variation  in  (10.15).  Therefore,  if  one  disre' ards  terms 

of  the  order  of  /m  w  /u>„  one  obtains  Eq.  (10.7).  One  may  reintro- 

m  M 

duce  the  time  dependence  of  the  power  into  the  linearized  form 
of  the  equation  for  frequency  (10.3  4).  Setting  1/[1  +  (P/P  ) ]  « 

o 

1  -  (P/Ps)  and  disregarding  the  product  5  P: 
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where 
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Integration  gives 
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where 


a 


The  effective  bandwidth  which  is  contained  in  the  pulse  formula 
(10.5)  starts  out  at  to,.  ,  the  value  prescribed  by  the  line 
width  in  the  absence  of  power,  and  approaches  asymptotically 
the  steady  state  value 
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Hence,  the  bandwidth  of  the  pulse  increases  as  a  function  of 
time  as  the  laser  medium  is  saturation  broadened.  The  bandwidth 
as  a  function  of  time  is  shown  in  Figure  9. 


XI.  Stability  of  Solutions  of  Homogeneously  Broadened  Laser  Medium 


In  the  subsequent  sections  we  study  the  stability  of  the 
supermodes  with  respect  to  arbitrary  perturbations.  It  is  as¬ 
sumed  that  the  mode  under  consideration,  say  the  n-th  super¬ 
mode,  is  perturbed,  the  perturbation  is  expanded  in  the  complete 
set  of  eigenfunctions  of  the  differential  equation  (harmonic 
oscillator  equation) .  Growth  of  any  one  component  of  the  per¬ 
turbation  signifies  instability. 

We  shall  follow  the  analysis  of  stability  of  Kurokawa. 

In  this  analysis,  the  perturbation  is  assumed  to  be  a  slow 
function  of  time  (compared  with  the  time  variation  of  the  un¬ 
perturbed  oscillation) .  This  condition  is  clearly  met  in  most 
laser  systems. 

In  the  same  way  as  we  treated  the  buildup  of  a  supermode, 
we  treat  the  time  evolution  of  a  perturbation  in  terms  of  its 
complex  frequency  6w ,  where  Im  6w  is  the  rate  of  decay  (if 
positive)  of  the  perturbation,  in  general  a  function  of  w  (i.e. 

and  hence  mode  number).  Because  we  assume  that  the  per¬ 
turbation  61 (w)  is  in  general  complex  and  time  dependent,  fre¬ 
quencies  of  oscillation  of  the  perturbed  modes  are  allowed  to 
deviate  from  the  steady  state  frequency.  In  other  words,  we 
treat  the  most  general  time  dependent  perturbation  with  the  only 
restriction  that  the  time  rate  of  change  of  the  perturbation  be 
slow  compared  to  the  cavity  mode  bandwidth.  Taking  a  perturba¬ 
tion  of  (10.2)  ,  and  taking  into  account  that  j 6 oj  x  1  61  = 
xc'  d/dt  61  and  assuming  a  resonant  modulation,  we  obtain 
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(11.1) 


Here  we  have  denoted  the  unperturbed  solution  by  the  subscript 
n,  assuming  that  the  n-th  supermode  exists  within  the  cavity. 
Next,  consider  the  perturbation  In6r. 


(11. 2  > 


If  the  current  I(w)  is  so  normalised  that  /|l(w)|2  dw  is 
equal  to  the  power ,  then 


6P  =  /lIn(b))  +  6I<WH2  "  / 1  Xn  (a>)  |  2  dw  «  2  /In(u)  Re  61  du>  (11.3) 
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It  is  convenient  to  express  61  as  a  superposition  of  the 
complete  set  of  eigenfunctions  of  the  harmonic  oscillator  equation: 
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The  eigenvalues  of  these  equations  are  (compare  4.9) 


E  =  2  /Hr  —  [n  +  - 


(11.5) 


The  eigenfunctions  are  identical  with  the  supermodes  (4.8).  Wc 

make  them  all  real  and  normalize  them  so  that  / u  u  dto  =  6  (11.6) 

n  m  nm 

We  get 
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and  introduce  this  expression  into  (11.1).  We  obtain 
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We  multiply  both  sides  by  um(w),  integrate  over  all  fre¬ 
quencies,  and  use  the  orthonomality  relation  (11.6).  The  result  is 
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dw  =  x  '  —  a 
c  dt  m 


(10.18) 


The  term  under  the  integral  couples  m  and  n  if  they  differ  by 
±1  or  ±  2.  Concentrating  on  m  =  n  -  1,  and  studying  only  the 
real  part  of  (11.8)  we  find 


(43) 


~,n  ,  -  E  )  Rc  a  .  =  x  '  —  Re  a  . 
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which  gives  exponential  arowth  because  E  ,  <  E  .  The  n-th 

n-l  n 

super-mode  is  unstable  with  respect  to  a  perturbation  proportional 
to  the  mode  of  the  next  lower  order. 

Next,  set  n  =  m  =  0  and  study  the  resulting  equation. 

One  has  one  equation  for  the  imaginary  part  of  a^,  another  one 
for  the  real  part 
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where  we  have  used  the  fact  that 
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The  component  Re  aQ  decays  at  the  rate 
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(11.12) 


(44) 


This  decay  rate  is  always  positive,  even  though  (11.12)  seems 
to  permit  negative  values.  We  recall,  however,  that  we  have  ex¬ 
panded  a  Lorentzian,  an  approximation  valia  only  when  the  pulse 
bandwidth  is  small  compared  with  If  the  expansion  had 

not  been  made,  (11.10)  would  be  positive  definite. 

The  component  Im  aQ  has  no  growth  or  decay.  Hence  a  cruad- 
rature  perturbation  aQ  uQ  (to)  experiences  no  growth  or  decay: 
i.e.  there  is  no  restoring  force  for  a  perturbation  of  this  kind. 

A  Quadrature  perturbation  corresponds  to  a  phase  perturbation  of 
the  carrier  frequency.  Hence,  the  mode  locked  pulses  have  no 

*  ■  V 

carrier  phase  stabilization,  just  like  a  free  running  van  der  Pol 
oscillator  [8] . 


(45) 


XII.  The  stability  of  Supermodes  for  Inhomogeneous  Laser  Medium 


Whereas  closed  form  solutions  do  not  exist  for  the  inhomo- 
geneously  broadened  medium,  it  is  still  possible  to  carry  out  a 
stability  analysis  of  the  supermodes  whose  aualitative  features 
have  been  determined  in  Section  IX.  The  differential  equation 
for  the  steady  state  solution  for  resonant  mode  locking  is 
(compare  9.1) 
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The  situation  in  the  present  case  is  simpler  than  for  the  homo¬ 
geneously  broadened  medium,  because  we  disregard  the  reactance 
associated  with  the  laser  medium.  We  write  down  the  equation 
governing  a  perturbation  61,  as  an  obvious  extension  of  the 

preceding  section.  Assuming  that  the  steady  state  is  the  n-th 
supermode : 
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We  may  separate  this  equation  into  real  and  imaginary  parts, 
obtaining  for  the  imaginary  part 
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Im  6l  =  —  Im  5l 
dt 


(12.3) 


This  equation  already  can  be  used  to  demonstrate  that  all 
higher  order  supermodes  are  unstable.  Indeed,  let  us  expand  the 
perturbation 


Im  61  *  E  a  u  (w)  (12.4) 

P  p  p 

where  the  up(w)  are  the  orthonormal  eigenfunctions  of  the 
differential  eauation 


(12.5) 


By  definition  =  -1  for  the  eigenfunction  corresponding  to 
the  steady  state  solution,  which  itself  is  an  eigenfunction  of 
(12.5).  Now  it  is  clear  from  the  nature  of  the  problem  that 
E  <  -  1  for  an  eigenfunction  with  less  peaks  (less  kinetic  en- 
eray  and  hence  less  negative  curvature)  than  the  assumed  unper¬ 
turbed  solution.  Introducing  (12.4)  into  (12.3),  multiplying  the 


equation  by  u^fa),  integrating  over  all  w  and  using  the  ortho- 
normality  condition,  one  obtains 
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(12.6) 


Equation  (12.6)  shows  that  a^  grows  exponentially  for  an 
eigenfunction  m  for  which  Em  <  "If  i.e.  for  an  eiqenfunction 
with  less  peaks  than  the  lowest  order  eigenfunction.  Hence,  all 
supermodes  are  unstable  except  (maybe)  the  one  with  a  single  peak. 
We  turn  to  the  analysis  of  this  mode,  which  we  denote  by  the  sub¬ 
script  0.  Consider  (12.6)  applied  to  a  perturbation  proportional 
to  the  lowest  order  eiqenmode,  aQ .  Since  E  *=  -  1,  we  find 

d 

—  Im  a.  =  0 
dt  0 


The  perturbation  experiences  no  "restoring  force".  This  is 
entirely  analogous  to  the  homogeneously  broadened  case  in  which 
quadrature  perturbations  also  did  not  experience  restoring  forces. 
Now  consider  the  eouation  for  the  real  part  of.  the  amplitude 
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We  may  expand  Re  6 1  in  terms  of  the  eiqenf unctions  un 
of  the  differential  eauation 
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This  set  of  eigenfunctions  does  not  contain  the  unperturbed 
steady  state  solution  because  of  the  additional  term  in  parentheses. 
It  is  clear  however  that  the  eigenfunction  of  (12.8)  with  a  single 
peak  denoted  uQ(to)  possesses  the  lower  i  eigenvalue  and  hence  is 
the  most  "dangerous"  so  far  as  growth  is  concerned. 

If  the  term  in  the  brackets  containing  Aq2 (w)  were  missing, 
u0(w)  would  reduce  to  the  steady  state  solution,  and  Eq  would 
be  equal  to  -1. 

It  is  clear  that 


by  itself  represents  the  potential  well  of  the  steady  state  problem  . 
The  additional  term  modifies  the  well.  The  deeper  the  well,  the 
lower  the  eigenvalue,  the  -higher  the  well,  the  higher  the  eigen- 
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value.  It  is  clear  that  the  additional  term  raises  the  well  and 
hence  raises  the  eigenvalue.  Therefore 

E  >  -1. 
o 

Returning  to  (12.8),  introducing 


Re  61  =  E  a  u  (w) 

P  p  p 

multiplying  both  sides  of  (12.8)  by  vQ  (w)  and  integrating  over 
all  a),  we  obtain 
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and  since  Eq  >  -1,  we  find  that  the  perturbation  aQ  uQ(u) 
decays.  Perturbations  Droportional  to  higher  order  eigenfunctions 
deca^  even  faster.  In  this  way  we  have  proven  the  stability  of 
the  lowest  order  supermode  of  the  inhomogeneously  broadened  case. 


XIII.  Conclusions 


The  replacement  of  the  discrete  cavity  spectrum  by  a  con¬ 
tinuum  made  the  treatment  of  mode  locking  in  the  frequency  domain 
a  relatively  simple  mathematical  problem  which  gives  physical 
intuition  a  free  reign.  We  were  able  to  find  all  the  steady  state 
solutions  of  the  homogeneously  broadened  laser  medium,  but  found 
that  only  the  lowest  order  solution  was  a  stable  one. 

The  inhomogeneous ly  broadened  laser  medium  does  not  have 
modes  which  can  be  evaluated  in  closed  form,  yet  their  physical 
properties  were  easily  deduced  from  the  particle  motion  analog. 
That  this  kind  of  approach  is  a  useful  one  was  confirmed  by  the 
stability  analysis  of  the  supermode  of  the  inhomogeneous ly  broad¬ 
ened  laser  for  the  execution  of  which  it  was  not  necessary  to 
find  the  detailed  functional  dependence  upon  frequency  of  the 
supermode.  Again  it  was  found  that  the  only  stable  sup-rmode  was 
the  one  with  a  single  peak  when  plotted  versus  frequency. 

The  effect  of  detuning  was  obtained.  The  build-up  of  the 
stable  supermode,  as  a  function  of  time,  was  evaluated  and  the 
attendant  change  in  bandwidth  of  the  transient  supermode  was  de¬ 
termined  for  the  case  of  a  homogeneously  broadened  laser  medium. 

It  is  believed  that  the  potential  of  the  present  analysis 
has  been  hardly  tapped  and  that  many  other  issues  of  interest, 
both  in  forced  mode  locking  and  saturable  absorber  mode  locking, 
will  be  analyzable  using  this  new  approach. 
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FIGURE  CAPTIONS 

1.  Equivalent  circuit  for  single  mode,  injection  locked 
oscillator. 

t 

2.  The  equivalent  circuit  for  multimode  injection  locked 
oscillator  system. 

3.  Eigensolutions  (Supermodes)  of  the  homogeneously  broadened 
laser  in  the  frequency  domain. 

4.  The  force-field-potential  as  a  function  of  "time". 

5.  The  normalized  amplitude  versus  normalized  frequency  of 
the  lowest  order  supermode  for  an  inhomogeneous ly 
broadened  laser. 
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6.  Range  of  normalized  gain  and  normalized  mode  locking 
modulation  within  which  steady  state  mode  locking  solutions 
with  single  peak  are  found. 

7.  Amplitude  at  line  center  as  function  of  normalized  modulation. 

8.  Transient  build-up  of  normalized  power  in  the  cavity. 

9.  Bandwidth  as  a  function  of  time. 
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